Mon. Not. R. Astron. Soc. 000, (1999) Printed 1 February 2008 (MN M^pC style file vl.4) 



Binary quasars 

Daniel J. Mortlock, 1 * Rachel L. Webster 1 * and Paul J. Francis 2 ' 3 * 

1 School of Physics, The University of Melbourne, Parkville, Victoria 3052, Australia 

2 Research School of Astronomy and Astrophysics, The Australian National University, Weston Creek, A.C.T. 2611, Australia 

3 Department of Physics and Theoretical Physics, Faculty of Science, The Australian National University, Canberra, A.C.T. 0200, 

Australia 



o 
o 
o 

< 



> 

On 
\D 
O 
OO 
O 
O 
O 

S3 

a,: 

6 

CO 

C3 



X 



Accepted 1999 June 15. Received 1999 June 15; in original form 1999 March 30 



ABSTRACT 

Quasar pairs are either physically distinct binary quasars or the result of gravitational 
lensing. The majority of known pairs are in fact lenses, with a few confirmed as 
binaries, leaving a population of objects that have not yet been successfully classified. 
Building on the arguments of Kochanek, Falco & Muhoz (1999a), it is shown that 
there are no objective reasons to reject the binary interpretation for most of these. In 
particular, the similarity of the spectra of the quasar pairs appears to be an artifact 
of the generic nature of quasar spectra. The two ambiguous pairs discovered as part 
of the Large Bright Quasar Survey (Q 1429—053 and Q 2153—0256) are analysed 
using principle components analysis, which shows that their spectral similarities are 
not greater than expected for a randomly chosen pair of quasars from the survey. The 
assumption of the binary hypothesis allows the dynamics, time-scales and separation 
distribution of binary quasars to be investigated and constrained. The most plausible 
model is that the quasars' activity is triggered by tidal interactions in a galatic merger, 
but that the (re-) activation of the galactic nuclei occurs quite late in the interaction, 
when the nuclei are within 80±30 kpc of each other. A simple dynamical friction model 
for the decaying orbits reproduces the observed distribution of projected separations, 
but the decay time inferred is comparable to a Hubble time. Hence it is predicted that 
binary quasars are only observable as such in the early stages of galactic collisions, 
after which the quiescent super-massive black holes orbit in the merger remnant for 
some time. 

Key words: quasars - galaxies: nuclei - galaxies: collisions - gravitational lensing. 



1 INTRODUCTION 

A small fraction (~ 1 in 1000) of quasars are observed to 
have at least one nearby quasar image with essentially the 
same redshift (e.g. Kochanek 1995; Hewett et al. 1998), and 
approximately 50 such systems are known in all (e.g. Keeton 
& Kochanek 1996) . Some of the pairs (and all of the higher 
multiples) are the result of gravitational lensing by interven- 
ing mass distributions; some are physically distinct binary 
quasars. In most cases the correct interpretation is quite 
clear, even if considerable observational effort was required, 
but there is a small number (~ 20) of pairs of quasars which 
cannot be categorised with certainty (e.g. Schneider 1993; 
Kochanek, Falco & Munoz 1999a; Peng et al. 1999). These 
tend to have angular separations of <; 3 arcsec, and in most 
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cases are treated as potential 'wide separation' lenses (e.g. 
Narayan & White 1988; Wambsganss et al. 1995; Mortlock, 
Webster & Hewett 1996; Park & Gott 1997; Williams 1997). 
If they are lenses they indicate the presence of a significant 
and otherwise unknown population of dark objects with the 
mass of groups or clusters of galaxies (e.g. Kochanek 1995; 
Hawkins 1997) . If they are binaries they are almost certainly 
the result of enhanced fuelling in galactic mergers or inter- 
actions (Djorgovski 1991; Kochanek et al. 1999a). 

The arguments in support of both possibilities are sum- 
marised in Sec. ^| In Sec. ^ principal components analysis 
is used to objectively test the spectral similarity of several 
quasar pairs. The binary hypothesis is subsequently adopted 
for the entire pair population, and the implications for the 
dynamics of binary quasars are explored in Sec. ^. The con- 
clusions are summarised in Sec. ^, together with a discussion 
of the future observational and theoretical possibilities. 
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Table 1. Possible binary quasars. 
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Status 
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MG 0023+171 


Binary? 


2 R 2 


1.2 





95 


170 + 150 


km 


s~ 


1 


4" 8 


29 


kpc 


34 kpc 


38 


kpc 


1 


Q 0101.8-3012 


Unknown 


O 2 


0.8 





89 


+ 200 


km 


s~ 


1 


n'.'o 


101 


kpc 


117 kpc 


132 


kpc 


2 


Q 0151+04480 


Binary? 


O 2 


3.6 


1 


91 


520 ± 160 


km 


s~ 


1 


3" 3 


19 


kpc 


25 kpc 


28 


kpc 


3 


QJ 0240-343 


Unknown 


o 2 


0.8 


1 


11 


250 ± 180 


km 


s~ 


1 


6'.'1 


37 


kpc 


46 kpc 


52 


kpc 


4 


Q 1120+01950 


Unknown 


o 2 


4.6 


1 


47 


200 ± 100 


km 


s~ 


1 


&.' 5 


40 


kpc 


49 kpc 


55 


kpc 


5 


PKS 1145-071 


Binary 


2 R 


0.8 


1 


45 


80 + 60 


km 


s~ 


1 


4'.' 2 


26 


kpc 


32 kpc 


35 


kpc 


6 


Q 1208+1011 


Lens? 


o 2 


1.5 


3 


80 


1200 ± 1200 


km 


s~ 


1 


CX'5 


2 


kpc 


3 kpc 


3 


kpc 


7 


HS 1216+5032 


Binary 


2 R 


1.8 


1 


15 


260 + 1000 

J 1 1UUU 






1 


9'.'1 


56 




69 kpc 


77 




8 


Q 1343+2640 


Binary 


2 R 


0.1 


2 


03 


200 ± 900 


km 


s~ 


1 


9'/5 


55 


kpc 


72 kpc 


79 


kpc 


9 


Q 1429-0053 


Unknown 


O 2 


3.1 


2 


08 


260 ± 300 


km 


s~ 


1 


5^'l 


30 


kpc 


39 kpc 


42 


kpc 


10 


Q 1634+267 


Unknown 


o 2 


1.6 


1 


96 


30 + 90 


km 


s~ 


1 


3"8 


22 


kpc 


29 kpc 


32 


kpc 


11 


J 1643+3156 


Binary 


2 R 


0.6 





59 


80 + 10 


km 


s~ 


1 


2" 3 


12 


kpc 


14 kpc 


16 


kpc 


12 


Q 2138-431 


Unknown 


O 2 


1.2 


1 


64 


+ 100 


km 


s~ 


1 


4" 5 


27 


kpc 


34 kpc 


38 


kpc 


13 


Q 2153-2056 


Binary? 


o 2 


3.4 


1 


85 


1100 ± 1500 


km 


s~ 


1 


7'.'8 


46 


kpc 


60 kpc 


66 


kpc 


14 


MGC 2214+3550 


Binary 


2 R 


0.5 





88 


150 ± 400 


km 


s~ 


1 


3'.'0 


18 


kpc 


21 kpc 


23 


kpc 


15 


Q 2345+007 


Lens? 


O 2 


1.5 


2 


15 


480 ± 500 


km 


s~ 


1 


7"3 


42 


kpc 


56 kpc 


61 


kpc 


16 



'Status' summarises the current evidence concerning the nature of each pair (See Sec. 2.1 for the classification criteria, and Mortlock 
(1999) for a summary of the evidence pertaining to each pair.); confirmed lenses are not included. 'Type' summarises the detected 
optical and radio emission of the pair using the notation of Kochanek et al. (1999a), where 2 R 2 denotes a pair in which both 
components are detected in the radio (and optical) ; 2 R denotes a pair with only one radio-loud component (which is hence a binary 
quasar); and O 2 denotes a pair with no radio emission at all. Am is the magnitude difference of each pair (in the optical); z is the 
redshift; |Ai>||| is the observed line-of-sight velocity difference; A8 is the angular separation of the two co mpon ents; and R is the 
projected physical separation of the two quasars, given for the three cosmological models described in Sec. [4.1.1 , with Ho = 70 km 
s — 1 Mpc -1 assumed throughout. 

References: 1. Hewitt et al. (1987); 2. Boyle et al. (1998); 3. Meylan et al. (1990); 4. Tinney (1995); 5. Meylan & Djorgovski (1988, 
1989); 6. Meylan et al. (1987) and Djorgovski et al. (1987); 7. Magain et al. (1992) and Maoz et al. (1992); 8. Hagen et al. (1996); 9. 
Crampton et al. (1988); 10. Hewett et al. (1989); 11. Djorgovski & Spinrad (1984); 12. Brotherton et al. (1999); 13. Hawkins et al. 
(1997); 14. Hewett et al. (1998); 15. Munoz et al. (1998); 16. Weedman et al. (1982). 



a Q 0151+0448 is also known as PHL 1222 and UM 144. 
b Q 1120+0195 is also known as UM 425. 



2 KNOWN QUASAR PAIRS 

Table [j] lists all known quasar pairs which have not been 
confirmed as gravitational lenses, in order of increasing right 
ascension. This list has been compiled mainly from previous 
lens candidate compilations (Turner 1988; Surdej 1990a,b; 
Surdej et al. 1991; Surdej 1991; Schneider, Ehlers & Falco 
1992; Surdej & Soucail 1993; Schneider 1993; Keeton & 
Kochanek 1996; Kochanek et al. 1999a), but also includes 



mor e re cent discoveries, such as Q 0101.8—3012 ( Boyle et al i 
1998^nd the small separation pair Q 1208+1011 (Magain 
et al. 1992; Maoz et al. 1992). The criteria by which these 
pairs are given their tentative classifications are discussed 
in Sec. 2.1, and a detaile d sum mary of the status of each 



pair is given in Mortlock ( 1999 ) . Some statistical arguments 



relating only to the overall properties of the population of 



pairs are presented in Sec. 2.2 



2.1 Classification criteria 

From the discovery of the first gravitational lens (Walsh, 
Carswell & Weymann 1979) and the subsequent interest in 
lensed quasars, it was soon realised that reasonably objective 
criteria must be established to determine whether pairs were 



lenses (e.g. Webster & Fitchett 1986; Schneider et al. 1992; 
Kochanek 1993; Kochanek et al. 1999a). 

A pair can only be positively confirmed as a binary if 
the spectra of the images are vastly different, if only one 
image is radio-loud (i.e. an '0 2 R' pair, in the notation of 
Kochanek et al. 1999a), or if the quasars' host are detected. 
The only necessary condition a binary must satisfy is that 
its components' spectra not be 'too similar' to eac h ot her, 
the meaning of which is investigated further in Sees. 2.2 and 



I 

The main sufficient conditions for a pair to be identi- 
fied as a lens are: the presence of more than two images^; 
the measurement of a time-delay between the images; or the 
detection of a plausible deflector. Naively, both a sufficient 
and necessary condition for a pair to be a lens is that its 
components' spectra are identical, as gravitational lensing 
is achromatic. However the difficulty in assessing just how 
similar two spectra are, combined with the generic nature 
of quasar spectra, means that it is difficult to prove a pair 
is not a binary on these grounds, as discussed further in 



t It is possible that a physical triple or quadrupl e qua sar exists, 
but this seems very unlikely, as discussed in Sec. [4.2.2|. 
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Sec. |§| Spectral similarity is not a necessary condition ei- 
ther - the light from a lensed quasar travels along different 
lines-of-sight to the observer, which can result in a num- 
ber of achromatic effects. Dust along one line-of-sight can 
redden individual images (e.g. Falco et al. 1999); microlens- 
ing by stellar-mass objects can magnify the continuum of 
one image relative to the emission lines; intrinsic variability 
coupled with the time-delay can result in the simultaneous 
spectra of components differing; and even the fact that the 
light that makes up each image comes from slightly differ- 
ent points in the quasar can result in different spectra being 
observed^} Further, the images of a lensed quasar can even 
yield different redshifts - the two components of Q 095 7+561 
have re dshifts which differ by Az = 0.0023 ± 0.00017 ( [Wills 
& WiU j 198C| ). 

Clearly the ambiguous quasar pairs do not satisfy any of 
these sufficient conditions, although some come close. The 
generic uncertain pair has no visible deflector, but quali- 
tatively similar spectra, neither of which allows a definite 
classification. 



2.2 The population of pairs 

If the 11 ambiguous pairs in Table [l] are predominantly 
lenses, the statistical properties of the sample should match 
those of the confirmed lenses, and likewise if they are mainly 
binaries. Their are some potential pit-falls to this approach 
(e.g. the definition of the samples; the presence of recently 
discovered pairs, such as Q 0101.8—3012, which may not re- 
quire the existence of a particularly dark lens) but it is a 
potentially powerful method . 

Kochanek et al. (1999a) used this technique to argue 
that most of the wide separation (i.e. A9 ^ 3 arcsec) pairs 
are not lenses for the following reasons. The existence of 
the 2 R binaries, combined with the knowledge that most 
quasars are radio-quiet, implies a comparable or larger pop- 
ulation of O 2 binaries - presumably the majority of the 10 
in Table [lj. The distribution of the ambiguous pairs' red- 
shifts is significantly different to that of the confirmed lenses, 
which is very difficult to account for in terms of lensing (e.g. 
Williams 1997). The distribution of flux ratios of the pairs 
and of the confirmed two- image lenses differ greatly (Fig. 0), 
although neither are consistent with the distribution pre- 
dicted by a simple, spherical lens model, implying that fur- 
ther investigation is warranted. Lastly, if the pairs are lenses, 
the absence of a smooth fall-off at larger separations is at 
odds with most theoretical models that predict wide sepa- 
ration lenses (e.g. Narayan & White 1988; Kochanek 1995). 

The statistical objections to the binary hypothesis have 
been the uncanny similarities of pairs' spectra (which is ad- 
dressed in Sec. |3J) and the high number of pairs at such 



t It is possible, for example, for a lens to have broad absorption 
features in only one component if the absorbing 'clouds' around a 
quasar are sufficiently small. The physical separation of photons 
in the two light paths would be only a few pc as they passed 
through the clouds, but it is possible that the broad line region 
consists of sh eets of plasma only a few metres thick (Blandford 
fe Rees 1991 ), in which case the two lines-of-sight are effectively 
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Figure 1. A cumulative plot of the magnitude differences, Am, 
between the possible lensed quasar pairs (dashed line), and all 
known two-image lenses, both with (upper solid line) and without 
(lower solid line) the known arcs and rings (which are generally 
treated as having Am = 0). The two smooth solid curves show 
the distribution predicted by a simple lensing model for survey 
magnitude limits of B = 19 (upper line) and B = 21 (lower line). 



small separations (relative to the predictions of the generic 
qua sar-qua sar correlation function). However, Kochanek et 
al. (1999a) showed that the number of pairs observed was 
consistent with the number expected from the observed 
small-scale quasar-galaxy correlation function. Further im- 
plications of this are explored in Sec. W. 



3 SPECTRAL SIMILARITIES OF QUASAR 
PAIRS 

The spectra of many of the pairs do appear unusually alike, 
but, before this can be regarded as an objection to the bi- 
nary hypothesis, it must be shown quantitatively that they 
are significantly more alike than those of randomly chosen 
quasars^} There have been several attempts to implement 
this sort of test, but the limitations of th e data have proved 
critical in most cases. Turner et al. (1988) found that the the 



unco rrclated. This possib ility is why O 0151+ 0448 (Mcylan et al 



1990) and Q 2153-2056 (Hcwctt et al. 19981 have been classified 



more tentatively here than by Peng et al. (199£). 



emission line shapes of Q 1634+267 A and B were more sim- 
ilar to each other than they were to those of another quasar 
of comparable luminosity and redshift, although the use of 
only a single compariso n ob j ect limited the strength of the 
conclusions. Peng et al. (1999) found that the C IV equivalent 
widths and continuum slopes of 14 pairs (all those in Table [j] 
bar Q 0101.8-3012 and Q 1208+1011) were more alike than 
over 97 per cent of comparable random quasar samples, im- 
plying tha t the apparent spectral similarities are real. Lastly, 
Hawkins (1997) found that the colours of the two compo- 
nents of Q 2138-431 and Q 2345+007 were closer than ex- 
pected, when compared to those of ~ 10 other quasars with 



S Even then, there are some reasons to think that binaries may 
have similar spectra due to th eir common environment and for- 
mation time ( Peng et al. 199s| ). However it is difficult to see how 
either could influence the localised velocity fields that are respon- 
sible for the emission line shapes. 
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similar redshifts. However the control quasars represent a 
heterogeneous sample, which can only increase the relative 
similarity of the pairs' colours. This illustrates why such 
relative tests are better at demonstrating pairs can be bi - 
nary objects. In a related investigation, Small et al. (1997) 
showed that the differences between the component spectra 
of Q 1634+267 and Q 2345+007 were no larger than the 
temporal variations in other quasars. Hence these pairs are 
consistent with being lenses, despite the differences in their 
spectra. However this represents only a necessary condition, 
and is not a sufficient condition for the lensing interpreta- 
tion. The choice of experiment is probably an artifact the 
greater interest in finding lenses. 

The analysis that follows is directed towards determin- 
ing whether the two ambiguous quasar pairs in the Large 
Bright Quasar Survey (LBQS) are c onsis tent with being bi- 
naries. The data is described in Sec. 



3.3 



3.1 



and the method of 

comparison discussed i n Se c. fS.2j. The results and interpre- 
tation are given in Sec. 



3.1 



Quasar pairs in the Large Bright Quasar 
Survey 



The LBQS is a sample of 1055 quasars with msj <; 18.5, 
selected spectroscopically and without any reference to mor- 
phology (Hewett, Foltz & Chaffee 1995). Due to both its size 
and its well-characterised selection effects, it is an excellent 
source of data for studying the similarities and differences 
between quasar spectra. 

As part of the survey, there was a systematic search for 
companions within ~ 10 arcsec of each quasar. The com- 
panion search has sufficient dynamic range (Am ~ 3.5) to 
easily pick out most secondary lensed images, and also to 
detect potentially dimmer binary companions. The search 
has yielded three pairs so far (listed in Table ^|) , and it 
is 'unlikely that further pai rs will be identified ' ( Hewett 
et al. 1998; ). Q 1009-0252 (jHewctt ct al. 1994) is a con- 



firmed lens; Q 1429—0053 ( Hewett et al. 198E ) is a reason - 
able lens candidate; and Q 2153—2056 ( Hewett et al. 199£ ) 
is probably a binary quasar. The spectra of all three pairs, 
taken concurrently in each case, are shown in Fig. H, and 
it is clear that there are significant differences between the 
primary and companion of both Q 1009—0252 (the lens) 
and Q 2153—2056. However, a more quantitative, and hence 
objective, means of testing the similarity of the spectra is 
needed to proceed further. 

3.2 Comparison of quasar spectra 

The most common method of analysing spectra (aside from 
visual inspection) is to measure the properties of individual 
features, such as emission lines and continuum properties. 
The weakness of this approach is that the classification of 
the features is usually subjective at some level. 

It is preferable to use a completely quantitative method, 
such as principal components analysis (PCA; Whitney 1983; 
Murtagh & Hecht 1987; Mittaz, Penston & Snijders 1990). 
PCA can be used on pre-determined features, but, for the 
reasons given above, it is more powerful to perform PCA 
on the raw spectra. This 'spectral' PCA has the advantages 
that it is a completely objective analysis, uses all the avail- 
able data, and is adept at dealing with low signal-to-noise 




X5 



Hi 



8000 



(b) c 


i 

j 1429- 


-0053 


- 
■ 




- 


i' u ''' 


i 


i 1 ! i t !,' , < f (Pi } 




4000 


5000 6000 7000 
X (A) 


- (c) C 


2153- 


2056 




!,/; 

"i ' 

h 








i 

! 1 1 . *i i i i ' 


! w",j 




4 y 1 

H 


'1 II 

i 





4000 



4500 



X 



5000 
(A) 



5500 



Figure 2. The spectra of the primaries (solid lines) and com- 
panions (dashed lines) of the three quasar pairs in the LBQS: 
Q 1009-0252 (a); Q 1429-0053 (b); and Q 2153-2056 (c). In 
each case, both spectra have been smoothed (with a smoothing 
length between 5 and 10 A) and normalised to have an average 
flux of unity, to aid in the comparison. Note that all spectra are 
in the observed frame, not the rest frame. 
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Table 2. Quasar pairs in the LBQS. 



Name 




Status 




m A 




m B 


Am 


z 


|A«||| 


A0 


Ref. 


Q 1009- 


-0252 


Lens 


V 


= 17.9 


V 


= 20.5 


2.6 


2.74 


150 ±200 kms" 1 


1'.'5 


1 


Q 1429- 


-0053 


Unknown 


R 


= 17.7 


R 


= 20.8 


3.1 


2.08 


260 ±300 km s" 1 


5'.'1 


2 


Q 2153- 


-2056 


Binary? 


B 


= 17.9 


B 


= 21.3 


3.4 


1.85 


1100 ± 1500 km s" 1 


7"8 


3 



'Status' summarises the current evidence concerning the nature of each pair (See Table |lj and Sec. 2.1); 
m\ and raj are the magnitudes of the primary and its companion, respectively, and Am is the difference 
between the two (The band is given in the table.); z is the redshift of the pair; |Au||| is the observed 
line-of-sight-velocity difference; and A8 is the angular separation of the two components. 
References: 1. Hewett et al. (1994); 2. Hewett et al. (1989); 3. Hewett et al. (1998). 



ratio spectra. With each spectrum treated as a vector in a 
multi-dimensional space, the sample of spectra represents 
a distribution in this space, the centroid of which is the 
mean spectrum. Subtracting the mean from each spectrum 
yields a set of points (centred on the origin) which represents 
the variation between spectra. The vector (or spectrum, of 
sorts) along which this distribution is most elongated is the 
first principle component (PC) of the sample. Subtracting 
it from all the already mean-subtracted spectra allows the 
procedure to be repeated, generating the subsequent compo- 
nents. The subtraction of only the first few PCs usually (i.e. 
if the PCA is successful) leaves a condensed, spherical distri- 
bution of points that is essentially a hyper-sphere of noise. 
Almost all the information is contained the in the first few 
components; t he d ata compression involved can be large. 
Francis et al. (1992) found that the LBQS spectra could be 
expressed in terms of coefficients of ~ 10 components, as 
opposed to several hundred wavelength bins. 



Following Francis et al. (1992), two subsets were taken 



from the LBQS for the PCA. Sample 1 contains 325 quasars 
(including all three pairs) with 1.4 < z < 2.7, and rest-frame 
spectra covering the range from 1400 A to 2200 A. Impor- 
tantly, this range covers several prominent emission lines, 
notably Civ (1549 A) and Aim/Cm] (1858 A/1909 A). 
Sample 2 has only 209 quasars (including the the lens and 
the wide separation pair Q 1429—0053), but, with 1.8 < 
z < 3.3, covers a sligtly bluer part of rest-frame spectrum: 
1180 A to 1780 A. It includes the C IV emission line and the 
Lya/Nv complex (1216 A/1240 A). The other wide sepa- 
ration pair, Q 2153—2056, could not be included in sample 
2, as the available spectra (Fig. |^ (c)) do not include the 
Lya/Nv blend. A number of quasars, including the stated 
pairs, appear in both samples. 

Spectral PCA was performed on both subsamples, gen- 
erating a set of PC spectra for each sample, and a set of 
coefficients for each quasar. As expected, most of the varia- 
tion between spectra was contained in the first few compo- 
nents, which are shown in Fig. |§] (for sample 1), and Fig. [| 
(for sample 2). T he com ponents are similar to those shown 
in Francis et al. (1992) (as expected, given the related set 
of input spectra), and their interpretation is also similar. 
The first three components of both samples are related to: 
the cores of the emission lines (PC 1; Figs. ^ and ^ (b)); 
the continuum slope (PC 2; Figs. ^| and ^| (c)); and broad 
absorption lines (PC 3; Figs. | and § (d)). The 4th PC of 
samples 1 and 2 relate to the wings of of the C IV and Ly a 
emission lines, respectively (Figs. || and ^ (e)). It is promis- 



ing that the C IV line is prominent in the first PC, as it is the 
most prominent difference between the spectra of several of 
the confirmed binaries (Q 12 16+5032 and Q 1343± 2640), 

and Peng et al. (Il999h found 



(1988) 



and both Turner et al. 
its properties revealing 

Whilst the PCA has enabled each spectrum to be spec- 
ified in terms of several numbers (i.e. the coefficients of the 
first few PCs), each spectrum is still a multi-dimensional 
vector. The data can be further compressed by defining a 
metric in the PC space. The PCA naturally scales the vari- 
ance in each component with its relative importance, so a 
simple Euclidean metric is chosen. The difference, D, be- 
tween the spectra of two quasars A and B is then given by 

1/2 



Dae 



Cb,h) 



(1) 



where Ca.u and Cb,ti are the coefficients of the nth principal 
component of spectra of A and B, respectively, and n max is 
the highest component used. The latter does not have an 
appreciable affect on the results, as long as the first few 
components are included, so n max = 4 was used. 

Two sets of D- values were created for each sample: the 
differences between all the spectra in the sample; and the 
differences between each of the pairs' primaries and the rest 
of the sample (which determines the fraction of spectra in 
the sample that are closer to the primary in question than 
its companion). The resultant cumulative plots are shown 
for sample 1 and sample 2 in Figs. |B| and ^j, respectively. 

3.3 Implications for the quasar pairs 

The two components of the confirmed lens, Q 1009—0252, 
are c lear ly different. This is not unexpected - as discussed in 
Sec. 2.1, there are a number of reasons why lensed spectra 



need not be the same. However, this analysis shows that 
they can actually be less similar than a two spectra chosen 
at random: Fig. |B| shows that 75 per cent of random pairs 
are more alike than the components of the lens. However the 
main difference between the two spectra of the Q 1009—0252 
is the reddening of the companion spectrum (Fig. |^ (a)), 
which is probably due to d ust along the line-of-sight to the 
fainter image (see Sec. 



2.1 



and Falco et al. 1999). 
The components of Q 2153—2056 are also more dis- 
similar than those of random pairs. Naively, the lesson of 
Q 1009-0252 implies that Q 2153-2056 cannot be rejected 
as a lens on this basis alone. The difference between the spec- 
tra of Q 2153—2056 cannot be so readily explained in terms 
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Figure 3. The mean spectrum (a) and the first four component 
'spectra' (b-e) of sample 1, plotted as a function of rest-frame 
(i.e. in the frame of the quasars) wavelength. The mean spectrum 
is normalised to an average flux of unity, and all the principal 
components have zero mean. The 3rd and 4th PCs have been 
inverted so that emission features are positive. The prominent 
emission lines are Civ (at ~ 1550 A) and the Alin/Cin] blend 
(at ~ 1900 A). 



Figure 4. The mean spectrum (a) and the first four component 
'spectra' (b-e) of sample 2, plotted as a function of rest-frame 
(i.e. in the frame of the quasars) wavelength. The mean spectrum 
is normalised to an average flux of unity, and all the principal 
components have zero mean. The 3rd and 4th PCs have been 
inverted so that emission features are positive. The prominent 
emission lines are the Lyo/N V blend (at ~ 1230 A) and Civ (at 
~ 1550 A). 



of lensing - indeed it appears that the companion is a broad 
absorption line quasar, which is only explained simply^] if 
they are two distinct objects. This distinction reveals one of 
the limitations of PCA - there is no physics or plausibility 
input, and so it cannot tell that some spectral differences 
(specifically Q 1009—0252) might be due to dust. In this 
instance, a continuum fit to the spectra would circumvent 
this problem. Due to both the spectral diff erence s and the 
absence of a visible deflector, Hewett et al. ( 1995 ) interpret 



Q 2153—2056 as a probable binary, a conclusion which is 
confirmed here. 

The results for Q 1429—0053 (which is a more proba- 
ble lens a priori) are more interesting. For the wavelength 
range covered by sample 1, the components of Q 1429—0053 
are more similar than ~ 85 per cent of random pairs. It 



Sec Sec. 2.1 for a complex explanation. 



also appears that Q 1429—0053 is a slightly unusual quasar 
spectrally, and so only ~ 9 per cent of the spectra are more 
similar to the primary than its companion. The sample 2 
results paint a slightly different picture. Q 1429—0053 B is 
no more similar to Q 1429—0053 A than a random quasar 
from the sample. This is despite the fact that Q 1429—0053 
A is at least slightly unusual - the dashed line and the 
solid line (representing the whole sample) are quite sepa- 
rated in Fig. ^. Such results are perfectly consistent with 
Q 1429-0053 being a lens (c.f. Q 1009-0252), but are 
also perfectly consistent with the two spectra having been 
drawn at random from the LBQS, and hence consistent 
with Q 1429—0053 being two distinct quasars. However, 
almost all the difference between the two components of 
Q 1429—0053 can be attributed to the relative strengths of 
the Ly a-Nv complex (appearing at ~ 4800 A in Fig. ti (b)), 
the equivalent widths differ by ~ 50 per cent. The possi bilit y 
of microlensing of the the secondary image (See Sec. 2.1.) 
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Figure 5. The cumulative distribution of th e d ifference between 
two spectra, D, shown for sample 1 (Sec. fc.2; ). The solid line 
is the distribution derived from all possible pairs of spectra in 
the sample; the other lines are from the differences between each 
of the primaries of the three pairs (Q 1429—0053: dashed line; 
Q 2153-2056: dot-dashed line; Q 1009-0252: dotted line) and 
the rest of the sample. The three vertical lines show difference 
between each primary and its respective companion. 
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Figure 6. The cumulative distribution of th e d ifference between 
two spectra, D, shown for sample 2 (Sec. fc.2; ). The solid line 
is the distribution derived from all possible pairs of spectra in 
the sample; the other lines are from the differences between each 
of the primaries of the two pairs (Q 1429—0053: dashed line; 
Q 1009—0252: dotted line) and the rest of the sample. The two 
vertical lines show difference between each primary and its re- 
spective companion. 



means that the lens interpretation is still possible, although 
the macrolensing flux ratio this would imply (Am ~ 3.5) is 
very unlikely a priori (e.g. Fig. [l]). 

There are several potential systematic errors which, 
whilst they do not affect the above conclusions, are im- 
portant in principle. The samples covered a wide redshift 
range, possibly increasing the spread in the samples' spec- 
tra; this could result in a true binary being rejected, but 
could not result in a lens being wrongly classed a binary. 



Another potential bias that is unimportant is the faintness 
of the companion images in these pairs - lower luminosity 
quasars might have different spectra. If the pairs are lenses 
this is clearly irrelevant, as they are not intrinsically so faint. 
If they are binaries, any bias that might make the spectra 
appear more dissimilar is acceptable, even desirable. 

The main source of random uncertainty is the low qual- 
ity of the companions' spectra (despite the fact that PCA is 
adept at analysing low signal-to-noise data) , which is simply 
a function of their faintness (Table ^). Several hours' ob- 
servation on a 10-m class telescope would yield companion 
spectra of better quality than the current primary spectra. 

On the basis of these results and the existing work sum- 
marised in Sec. |i[ all the ambiguous pairs are assumed to 
be binaries in Sec. ^. This, however, is quite an adventurous 
extrap olatio n - the broader but less rigorous results of Peng 
et al. imply that the spectra of the pair population 

are unusually alike. 



4 DYNAMICS OF BINARY QUASARS 

If all the quasar pairs listed in Table [l] are physical bina- 
ries, there are far too many to represent a simple a smal l- 
scale manifestation of quasar clustering (Djorgovski 1991). 
The fairly natural conclusion (Djorgovski 1991; Schneider 
1993) is that any close quasar pairs are the result of in- 
cre ased nu clear activity in interacting galaxies. Kochanek et 
al. (1999a) showed that the number of binaries matches the 
number expected from the quasar-galaxy correlation func- 
tion at small scales. Henceforth, all the pairs in Table |l| are 
treated as physical binaries, and the merger model is also 
adopted. The sample of pairs can then be use d to constrain 
the distribution of physical separations (Sec. |4.l| ), and also 
to investigate the formation and evolution of binary quasars 
(Sec. O). 



4.1 Physical separations 

The projected separations of the qu asar p airs can be inferred 
for a given cosmological model (Sec. 4.1.1 ) , but some method 
of inversion is required to obtain the distribution of ph ysical 



separations. A random deprojection is explored in Se c. 4.1.2 
and a simple orbital model is discussed in Sec. 4.1.3. 



4-1.1 Cosmological models 

The projected separations of the pairs are given by R = 
A8 oIa(0, z), where d\ is the angular diameter distance, and 
so R varies with the cosmological model (and also within 
a model, due to weak lensing). Three models are used: the 
Einstein-de Sitter (EdS) model (f2 mo = 1 and £Ia = 0); a 
low-density model (f2 mo = 0.3 and Q\ = 0); and a low- 
density, spatially-flat model (fim = 0.3 and flA = 0.7), 
where fi mo is the ratio of the current density of the universe 
to the critical density and Qa is the similarly normalised 
cosmological constant (e.g. Carroll, Press & Turner 1992). 
The ~ 10 per cent error in Hubble's constant (taken to be 
Ho — 70 km s -1 Mpc -1 here) does not contribute signifi- 
cantly to the uncertainty in R. 
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Figure 7. The di fferen tial (a) and cumulative (b) separation distribution of binary quasars obtained by applying the random deprojection 
describ ed in Sec. [4.1,2 to the data listed in Table l| The three lines show the results assuming the three cosmological models described 
in Sec. [4.1, l[ the EdS model (solid line); the low censity model (dashed line) and the low density, flat model (dot-dashed line). 



4-1.2 Random deprojection 

The problem of deprojection - effectively obtaining three 
dimensional information from a two-dimensional image of 
the sky - is an old one, with a number of different par- 
tial solutions (e.g. Courant & Hilbert 1962; Gerhard & Bin- 
ney 1996). The problem at hand is to obtain a probability 
distribution of the physical separation, dp/dr\n, given an 
observed projected separation, R. A Bayesian approach is 
adopted whereby the separation vector, r, is considered to 
be randomly-oriented, a posteriori. Mathematically, dp/di = 
sin(i), where i is the (inclination) angle between r and the 
line-of-sight. It is normalised for the range < i < it/2, and 
the projected separation is given by R — r|sin(i)|. From 
this, a change of variables yields 



dp 




dp 




di 


dr 


R 


d? 




dr 



16 R J 
tv r 4 



if r <R, 
if r > R, 



(2) 



where the final distribution is also normalised to unity, and 
has an expectation value of 16/(37r)_R, assuming R > 0. 

Convolving dp/dr|_R with the data gives an approxima- 
tion to the physical separation distribution of the binary 
quasars, the differential and cumulative forms of which in 
Fig. [t]. This 'random deprojection' shows a tail of pairs with 
quite large separations and a possible paucity of small sep- 
aration pairs. (The one pair with r & 10 kpc, Q 1208+1011 
is a probable lens; if it were confirmed as such, the small 
separation 'hole' would be quite clear.) Any hole cannot 
explained by selection effects, as there are a large number 
of lenses with angular separations of ~ 1 arcsec, implying 
that a ny suc h binaries would have been found. Kochanek 
et al 



(1999a 



showed that dynamical friction (Sec. i.2.c) 
can account for the hole, although the inferred separation 
distribution (linearly rising with r) is otherwise quite dif- 
ferent from the approximately exponential fall-off apparent 
in Fig. 0. Some clues as to the nature of small separation 
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Figure 8. The cumulative distribution of projected separation, 
R, for the quasar pairs in T able h i and the best-fit dynamical fric- 
tion model described in Sec. 4.1.3[ The three lines show the r esults 
assuming the three cosmological models described in Sec. 
the EdS model (solid line); the low density model (dashed line) 
and the low density, flat model (dot-dashed line). 



binaries might be given by the discove ry of three possibl e 
binaries in a sample of ~ 100 BL Lacs (Scarpa et al. 1999). 



The frequency of these pairs is far too high to be explained 
in terms of gravitational lensing by galaxies, but the binary 
interpretation is no more comfortable. It would require not 
only that a few per cent of BL Lacs are binaries, but that 
the optical jets of the companion BL Lacs are parallel. 

4-1.3 Orbital models 

An approach alm ost opposite to the random deprojection 
used in Sec. 1.1.2] is to assume an orbital model for the bi- 
naries, and fit the resultant distribution of projected sepa- 
rations to the data. The paths of two quasars in a pair of 
merging galaxies is undoubtedly complex, but a natural first 
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Figure 9. A scatter plot showing the projected separation, R, 
and the line-of-sight velocity difference, Aui | , of the candidate bi- 
nary quasars listed in Table The velocity error bars are random 
measurement errors, but the separation errors are taken from the 
spread in R over cosmological models, and are thus both random 
and systematic. Also shown is -B m i n = [as defined in equa- 
tion 
line). 



) , if each quasar's host galaxy weighed lu Mq (solid 



IO^Mq (daohcd line) or lO^M^ (dot daahcd line). If 

the points should lie below the relevant 



the binaries are bound, 
line. 



approximation is to assume that they are in decaying ellip- 
tical orbits. However averaging over the random orientation 
removes most of the effe cts of ellipticity , so circular orbits 
can be used in general (Mortlock 1999). The main differ- 
ence was in the radius of closest approach, wh ich m ay be 
relevant to the activation of the quasars (Sec. |4.2.2| ). The 
existence of the outer cut-off of between 50 kpc and 100 kpc 
is confirmed, independent of the distribution of separations 
used. The putative inner cut-off is less well constrained, and, 
assuming dp/dr oc r (based on the dynamical friction ap- 
proximation discussed in Sec. 4.2.3), the data is consistent 
with a population of pairs extending to zero separation, as 
shown in Fig. |§| However, formally better fits are obtained 
with an inner cut-off of ~ 10 kpc, especially if the smallest 
separation pair, Q 1208+1011, is a lens as suspected. An in- 
ner cut-off would imply that the at least one of the quasars 
in each pair 'turns off' prior to the complete decay of th e 
orbits. Possible reasons for this are discussed in Sec. i. 2. J . 



4.2 Physical processes 

If binary quasars are the result of galactic collisions during 
which both nuclei become active, the assumption that the 
syste ms are bound implies a minimum mass for the systems 
(Sec. 4.2.1), and the distribution of physical separations 
(along with timing arguments) can be used to invest igate 
bot h the activation and evolution of the pairs (Sees. 4.2.2 



and 4.2.3 



respectively). 



4-2.1 Total energy 

It is possible that quasars could form in 'glancing' colli- 
sions (e.g. Noguchi 1987, 1988), but the absence of 'post- 
collisional' pairs with larger separations argues against this 



for the binary population, implying that they exist in gravi- 
tationally bound systems. The minimum possible centre-of- 
mass frame energy of a pair (obtained by setting the line- 
of-sight separation and projected velocity difference to zero) 



M A M B 



A«j| GM A M B 



M A + M B 2 



R 



(3) 



where G is Newton's gravitational constant, and Ma and 
Mb are the masses associated with the two quasars. Invert- 
ing equation (J3J) and assuming Em'm — gives the hard lower 
limit for the mass of a bound system as 



Mtot = M A + Mb > 



RAvf { 
2G 



(4) 



Fig. [] shows a scatter plot of the projected separations 
and velocity differences of all the pairs in Table [l], with the 
large uncertainties in both R and Awn clearly important. 
Despite this, Mtot iS 5 x 10 11 Mq is inferred, confirming 
that quasars are typically in larger galaxies (Hooper, Impey 
& Foltz 199 7). Ho wever using the random deprojection dis- 
ussed in Sec. 4.1.2 to account for the projection effe cts yields 



a lar ger (but softer) limit of M tot ^ 2 x 10 Mq ( pVIortlock 



199S) 



4-2.2 (Re-) activation 

It is widely believed that quasars are the active nuclei of 
galaxies, where gas (and other material) accretes onto a mas- 
sive central black hole (e.g. Lynden-Bell 1969; Rees 1984). 
Many galaxies have massive but quiescent black holes at 
their centres (e.g. Kormendy et al. 1996, 1997), and there 
are several plausible explanations for how they may begin 
accreting and become active galactic nuclei. One impor- 
tant mechanism appears to be encounters between galaxies, 
which can result in the formation of a massive black hole 
(e.g. Rees 1984) or provide fuel to previously quiescent nu- 
clei. This is supported by Hubble Space Telescope (HST) ob- 
servations which show that many quasars' host galaxies are 
distorted (Bahcall, Kirhakos & Schneider 1995a; Boyce, Dis- 
ney & Bleaken 1999) and that low-redshift quasars are often 
found to be physically associated with other galaxies (e.g. 
Yee 1987; Bahcall, Kirhakos & Sc hneider 1995b; Bahcall e t 
al. 1997) or involved in collisions (Steidel & Sargent 1990). 
TV-body simulations also show that gravitational torques be- 
tween colliding galaxies can cause gas to flow to thei r cores, 
potentially acting as fuel for the central black holes (Barnee 
&i Hernquist 1996). The scales of the in-flows are not great 



enough to produce new black holes, but are consistent with 
the re-fuelling model. However, the simulations only predict 
in-flow to within a few hundred pc of the galaxies' cores, 
whereas t ypical quasar accretion disks extend out to only 
pc scales (Peterson 1997). It is not clear whether dissipation 



is sufficient to allow the gas to fall further. 

The formation of binary quasars is qualitatively consis- 
tent with the reactivation process, provided that the progen- 
itors (or their central black holes) are sufficiently massive. 
The observational result that the ac tivation radius is be- 
tween 50 kpc and 100 kpc (Sec. |4.l| ) implies that quasar 
formation does not occur until quite late in the collision 
process. This is consistent with the TV-body simulations of 
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Barnes & Hernquist (1996), which showed that gas inflows 



took some time to occur. 

The separation scale ol the binaries, together with the 
requirement ol a massive black hole, is also suggestive ol the 
lact that quasars form in larger galaxies. II this is the case, 
encounters between dwarfs and large galaxies can result in 
the formation ol a quasar, but only collisions between two 
large galaxies can result in a binary quasar. The traction ol 
merger-related quasars which are binaries is comparable to 
the fraction ol collisions in which both galaxies hav e cen - 
tral black holes ol <; 10 6 M Q . Assuming a Schechter ( 197Q 
mass lunction for the galaxies, and that ~ 20 per cent ol 
t hose heavier than M* ^ 10 Mm contain large black holes 
( Magorrian et al. 1998 ), about 1 per cent ol all merger- 
formed quasars should be in binary systemsjjj]. Combined 
with the observation that only 1 quasar in ~ 1000 is a bi- 
nary, this implies that about 10 per cent ol quasar activity 
is related to galactic interactions. This is probably an un- 
derestimate ol the fraction, however, as interactions between 
two large galaxies are almost certainly more efficient than 
collisions involving a dwarf galaxy. Further, il quasars all 
have a dusty torus which is both optically thick and sub- 
tends a large fraction ol the solid angle around the central 
engine, both the quasars in a pair must be close to lace-on 
for it to be observed as a binary. Extrapolating from the 
observation tha t ~ 25 per cent ol all Seylerts are ol Type I 
(Peterson 1997), 1 in ~ 4 ol all quasars are visible, but only 



1 in ~ 16 binaries are detected as such. It is thus possible 
that ~ 50 per cent ol all quasars are formed during galactic 
interactions and collisions. 



4.2.3 Evolution 

The theor etical and observational evidence discussed in 
Sec. 1.2.2] is consistent with binary quasars forming dur- 
ing galactic mergers. Once the merger has begun to set- 
tle, the quasars can be treated as point-masses moving in 
a static, isothermal halo. Their orbits slowly decay, as they 
lose energy and angular momentum through dynamical fric- 
tion (e.g. Binney & Tremaine 1987), eventually either merg- 
ing (e.g. Makino & Ebisuzaki 1994) or forming a hard binary 
that lasts for a Hubble time or more (e.g. Begelman, Bland- 
ford & Rees 1980; Rajagopal & Romani 1996). II the quasars 
have mass M ~ 10 9 Mq by this stage, and are in circular or- 
bits ol speed \/2cr (where a is the velocity dispersio n ol the 
halo) , the time sp ent with separation at radius r is (Binney 
fe Tremaine 1987| ) 

&i 3.3 _ r , (5) 



dr 



GM ln(A) 



where ln(A) ~ ln(2r*a 2 /GM) is the Coulomb logarithm, 
which charact erises the strength ol the interaction, and 
r* — 10 kpc (Lacey fc Cole 1993). Hence Ap/Ar oc r, as 



used in Sec. 4.1.3 to fit the observed separation distribution. 



The orbital decay proceeds at a greater rate as r decr ease s, 
potentially explaining the small separation hole (Sec. |4.l|). 



II It also implies that only ~ 1 per cent of multiple quasars should 
be triples etc., justifying the assumption that all observed triples 
and quadruples are actually lenses. 



The time-scale for the decay, £df, is found by inte- 
grating equation (|^) from r max to 0. Following Binney & 
Tremaine ( |l987[ ), 



toF 



1.65 ri 



ln(A) GM 



4 x 10 J 



1O 9 M 



20 kpc I 200 km s" 



M 



yr, (6) 



where the scale values have been chosen to give the short- 
est plausible decay timej^j. The dynamical time-scale ol the 
binaries is then close to a Hubble time, and much longer 
than both the 'settling' time ol the merged halo (< 10 9 yr; 
Barnes 1992) and the expected quasar lifetime. A black hole 
of initial mass ~ 10 7 A/q (as implied by dynamical measure- 
ments ol local galaxi es; Kormendy et al. 1996) accreting at 
the Eddington limit ( |Peterson 1997| ) would reach ~ 10 10 M Q 
in much less than 10 M yr. No such massive black holes are 
observed now (e.g. Kormendy & Richstone 1995), so this 
places an upper limit on quasar masses and hence their life- 
times. Comparison ol the three time-scales implies that bi- 
nary quasars are short-lived, existing only whilst the host 
galaxies are actively merging. The most likely explanation 
for the death ol the binaries (and other quasars in merg- 
ers) is that the in-flow ol gas ceases, probably as the merger 
becomes more stable. Moreover, by the time the orbits de- 
cay due to dynamical friction, the merging black holes are 
already quiescent, and the observed binary separation dis- 
tribution is a random sampling ol the chaotic phase ol the 
merger. 

II the above model is correct, the hosts ol binary 
quasars should be distinct (in the case ol a 'new-born' bi- 
nary), or, more likely, highl y distorted. The h ost galaxies ol 
MGC 2214+3350 A and B ([fvlunoz et al. 1998[) have been de- 



tected by the HST, and in lact o ut-shine the quasars in the 
//-band (Kochanek et al. 1999b). The galaxies show no ob- 
vious signs ol being distorted, which implies both that their 
physical separation is considerably greater than their pro- 
jected separation ol ~ 20 kpc, and that the nuclei have only 
recently become active. Comparable observations - which 
need not be prohibitively deep - ol some ol the other am- 
biguous pairs ought to prove similarly revealing. 



5 CONCLUSIONS 



Djorgovski (1991) and Kochanek et al. (1999a) have both 



argued that most of the known wide separation quasars pairs 
are not the result of gravitational lensing, as summarised in 
Sec. ^| However, there have also been some valid objections 
to the interpretation of the pairs as physical binaries. The 
strongest was that the spectra of the pairs appear to be 
very similar. Sec. |^ consists of a quantitative evaluation of 
this claim for the three quasar pairs discovered as part of 
the LBQS, and shows that, despite appearances, none of 



** The val ue of r max used is half the activation radius inferred 
in Sec. |-t. l| . as it is possible that both quasars orbit the centre 
of the halo, whence their distance between them is twice their 
orbital radius. Also, the black holes themselves might start with 
Al <C 10 9 Mq, but they must be associated with most of their 
eventual mass; hence the choice of the canonical mass. 
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these pairs can be rejected as binaries on spectral grounds. 
Extrapolating to the rest of the ambiguous quasar pairs, the 
simplest conclusion is that they are all binary quasars. 

The adoption of the binary hypothesis for the entire 
population of 16 pairs listed in Table l] provides a significant 
data-set of binary quasars, and even if some do turn out to 
be lenses, the statistical properties of the population will not 
be greatly changed. Assuming quasars are the nuclei of large 
galaxies, the velocity differences of the pairs are consistent 
with them belonging to bound systems, which supports the 
idea that they are part of galactic mergers. Bot h the orbital 
and deprojection analyses presented in Sec. 4.1 suggest that 



the galactic nuclei become active at separations of between 
50 kpc and 100 kpc, depending on the cosmology and the 
ellipticity of the quasars' orbits. This is certainly consistent 
with the scales at which the tides from large galaxies become 
important. However the time-scales suggested for the orbits 
of the nuclei to decay are very large - certainly longer than 
the lifetime of a quasar accreting at the Eddington limit. 
Indeed it seems probable that the quasars remain active only 
whilst the host galaxies are merging, in which case most of 
the binaries should exist in very distorted hosts, and not 
in a relaxed merger. The obser vation that the host gal axies 
of MGC 2214+3550 A and B ( [Kochanek et al. 1999bj ) are 
undisturbed then suggests that this is a very young binary. 

Whilst it is clear that the physics of merging galaxies 
and active galactic nuclei are far from solved, there are also 
a number of questions that can be answered observationally. 
Further spectroscopy and radio observations of a number of 
the quasar pairs should reveal whether they are binaries^ 
although there are also several which are still not classified 
despite major observational efforts. The nature of the 'turn 
on' radius will be mo re strongly constra ined by the 2 degree 
Field quasar survey (Boyle et al. 1998), which will include 



analysis of the nearby companions of ~ 3 x 10 qu asars. On 



a longe r time-s cale, the Sloan Digital Sky Survey ( Gunn & 
Weinberg 1995) will observe ~ 10 5 quasars, and is expected 
to find several hundred lenses and pairs, an order of mag- 
nitude increase on the number known presently. Both the 
lower limits on mass of the quasars' host galaxies and the 
constraints on the activation radius could be determined to 
within ~ 10 per cent from such a sample. 
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